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Abstract 

We study neutralino production at the linear collider with the subsequent two- 
body decays Xn^^ ^^^d —>■ ££, with £ = e,fi, r, or Z° — > qq with q = c,b. We 
show that transverse electron and positron beam polarizations allow the definition 
of unique CP observables. These are azimuthal asymmetries in the distributions 
of the final leptons or quarks. We calculate these CP asymmetries and the cross 
sections in the Minimal Supersymmetric Standard Model with complex higgsino 
and gaugino parameters fj, and Mi. For final quark pairs, we find CP asymmetries 
as large as 30%. We discuss the significances for observing the CP asymmetries at 
the International Linear Collider (ILC). Finally we compare the CP asymmetries 
with those asymmetries which require unpolarized and/or longitudinally polarized 
beams only. 



1 Introduction 



Supcrsymmctric (SUSY) models predict new particles with masses of the order of a few 
hundred GeV [1,2]. Their discovery is a major goal of present and future colliders in 
the TeV range. In particular, the International e+e^ Linear Collider (ILC) [3-7], with a 
center-of-mass energy of ^/s = 500 GeV and an integrated luminosity of C = 500 fb^^ in 
the first stage, will precisely measure the masses and couphngs of the SUSY particles when 
they are kincmatically accessible. It has been shown that the underlying parameters of 
the SUSY model can be determined at the percent level and better [3,7-9]. In particular, 
the option of polarized e'^ and beams [10] at the ILC can yield higher statistics to 
test models beyond the Standard Model (SM). Transversely polarized beams allow to 
study additional observables which are sensitive to effects of new physics. These are, for 
example, models with extra spacial dimensions, specific triple-gauge boson couplings, and 
also new sources of CP violation [10]. 

In the Minimal Supersymmetric Standard Model (MSSM) [1,2], the spin-half super- 
partners of the neutral gauge and CP-even Higgs bosons mix and form the four neutralinos 
At tree-level, the neutralino sector of the MSSM is defined by the U{1)y and SU{2)l 
gaugino mass parameters Mi and M2, respectively, the higgsino mass parameter /x, and 
the ratio tan/3 — vijvx of the vacuum expectation values of the two neutral Higgs fields. 
Besides the sleptons, the supcrpartncrs of the leptons, and the charginos, the superpart- 
ners of the charged gauge and Higgs bosons, the neutralinos are expected to be among the 
lightest SUSY particles in many models. The neutralinos will be pair-produced [11-13] 
at the ILC 

e+e- ^ x-X?, i,i = l,...,4. (1) 

By measurements of the neutralino masses, cross sections and decay distributions, meth- 
ods have been developed to determine the parameters of the neutralino sector [14-19]. If 
CP is violated, the parameters Mi = |Mi| e^^^^^ and \x = e*'^'" can be complex, while 
M2 and tan j3 can then be chosen real and positive. 

In general large values for CP phases in SUSY models lead to theoretical predictions 
for the electric dipole moments (EDM) of electron, neutron and that of the atoms ^^^Hg 
and 205rp|^ which are close or beyond the current experimental upper bounds [20-22]. 
The restrictions on the phases from EDM measurements, however, strongly depend on 
the SUSY model, sec e.g. [23,24], and on the scenario [25,26]. This means that SUSY 
CP phases of the order one are not ruled out by the present EDM experiments. For an 
unambiguous determination of the CP phases independent measurements are necessary, 
e.g., by analyzing CP sensitive observables at coUiders, in particular at the ILC. 

CP asymmetries with triple products [27, 28] have been analyzed for the production 
of neutrahnos and for their various two-body [29-34] and three-body decays [18,35-37]. 
It has been pointed out that longitudinally polarized beams can enhance simultaneously 
the cross sections and the CP asymmetries, such that higher statistics allows to measure 
even small phases [10]. 
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In contrast to longitudinally polarized beams, the possibility of transverse beam po- 
larization allows us to define a whole class of new and unique CP sensitive observables. 
These are asymmetries in the azimuthal distributions of final state particles. In general 
the transverse polarization of both beams is required, since these asymmetries vanish if 
only one beam is polarized. CP observables for transversely polarized beams have been 
studied for the production of selectrons [38] and charginos [39,40]. For neutralino pro- 
duction, two-body decays have been analyzed: x° '^lr'^^ ~^ ^~^^~Xi = [34,41], 
X°-X?/^°andx°-x;^M41]. 

In this work we study CP asymmetries with triple products in neutralino produc- 
tion (1), followed by the decay of one of the neutralinos into the Z° boson 

X°^Xn^°, n=l and j = 2, 3 , (2) 

for longitudinal and transverse beam polarizations. Due to the Majorana properties of the 
neutralinos, the angular distribution of the boson is independent of the polarization, 
if the polarization is summed [29,33,41-43]. Thus for longitudinally polarized beams, 
all CP sensitive information is lost. On the other hand, the CP sensitive observables in the 
reaction (1) with transversely polarized beams require the reconstruction of the neutralino 
production plane. This can be challenging, since there is not enough kincmatical infor- 
mation from energy and momentum conservation in the process e+e^ XiXi^^ [34,41]. 
Hence in either way, the polarization of the Z^-' has to be included and analyzed by the 
angular distributions of its decay products. On account of this, we study its leptonic and 
hadronic decays 

^°^//, / = ^,9, •^ = e,/i,r, g = c,6. (3) 

We will show that spin correlations between production and decay lead to sizeable CP 
asymmetries, which do not require the reconstruction of the production plane. 

In Section 2, we give the Lagrangians and couplings for neutralino production (1) and 
decay (2), (3). In Section 3, we present the analytical formulae for the amplitude squared 
with polarized beams. We define the CP asymmetries for transversely and longitudinally 
polarized beams in Section 4. In Section 5, we present numerical results for the CP 
asymmetries and the cross sections. We compare the CP asymmetries with longitudinal 
and transverse beam polarizations, and show that transverse beam polarizations can help 
to determine the phases in the neutralino sector. In Section 6, we give a summary and 
conclusions. 
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2 Lagrangians and couplings 

In the MSSM, neutralino production (1) proceeds via Z° boson exchange in the s-channel, 
and selectron eL,R exchange in the t- and u-channels, see the Feynman diagrams in Fig. 1. 



The Lagrangians for neutrahno production e^e — > x°x5 and decay x 
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with Pl^r — (1 ^75)/2. In the photino, zino, Higgsino basis (7, Z, H^, H^), the couplings 
are 

1 



2 



{NisN*, - NuN*,) cos 2/3 + {NisN*, + NuN*,) sin 2/3 



1 



- sin^ 6'vF)^i2 + sin OwNn 



.cos ^2 
f^i = V2 sin [ tan OwN*^ - N*^ 
Lf = T^f — Qf sin^ 9w , Rf = —Qf sin^ 9w , 



(7) 
(8) 
(9) 
(10) 
(11) 



with the weak mixing angle 6w, the weak coupling constant g = e/ sin^^y, e > 0, the 
electric charge q/ and isospin T^f of fermion /, and the ratio tan/3 = V2/V1 of the vac- 

"L R 

uum expectation values of the two Higgs fields. The neutralino couplings O^j ' and 
/^'^ contain the complex mixing elements Nij, which diagonalize the neutralino matrix 
^ioXafi-^pk ~ ""^xi^ik [1]; with the ueutralino masses m,^. > 0. 
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3 Distributions and cross section 



The CP sensitive observables for transversal and longitudinal beam polarizations are 
asymmetries in the angular distributions of final state particles. They depend on the spin 
correlations between neutralino production and decay (l)-(3). The polarizations of the 
particles can be included using the spin density matrix formalism [44]. In this formal- 
ism, the amplitude squared for neutralino production has been calculated at tree level 
for unpolarized and longitudinally polarized beams [12], and for transversely polarized 
beams [34]. The subsequent decay of one of the neutralinos into a Z° boson, x° X^^^i 
followed by fj\ has also been calculated in the density matrix formalism [29,33,43]. 

For completeness, we shortly summarize the results in the following. 

The amplitude squared for neutralino production (1) and decay (2)-(3) is given by 



\T\' = A\Am\'\A(Z' 



0m2 



6=1 



(12) 



with the propagators of the decaying neutralino and boson 



= - -T^T-r^^ A(z°) = - -3^-——, (13) 



where Fz and F^. are the total decay widths of the associated particles and sz and s-^^. 
are their invariant masses. The differential cross section then reads 

dcr = ^\T\'^dUps{s;p^.,p^„,pf,Pf) , (14) 

where liLips is the Lorcntz invariant phase-space clement, given in Appendix B. 

The amplitude squared (12) has contributions from neutralino production (P) and de- 
cay (D). The function P is independent of the polarizations of the neutralinos, whereas 
Ep depends on the polarization of neutralino x°. The longitudinal polarization of the 
neutralino is given by Sp/P, the transverse polarization of the neutralino in the produc- 
tion plane is given by Sp/P and the polarization perpendicular to the production plane 
is given by Sp/P. It is convenient to decompose the functions P and Sp into unpolarized 
(unpol), longitudinally (L) and transversely polarized (T) contributions of the electron 
and positron beams: 

P = Punpoi + Pl + Pt , Sp = Sp^j^pQi + T^p^^ + Sp,^ . (15) 

The quantities Punpoh Pl and Sp^j^p^jj, Ep^;^ are given in [12], and the quantities Pt and 
Spj, have been calculated in [34]. 

The quantities Sp^. with transversal beam polarization depend on the polarization of 
the neutralino Xj and are^ 

E^^ = E^(Zei)T + EUZg«)T + E^(eiefl)T , (16) 



'I'These terms are also given in [34], however, there in the terms ^p{ZeL,R)T and P{ZeL,R)T a factor 
g"^ / cos^ 9w is missing. 
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where 

Z COS uw 



Im(A(Z)/,t40:f [A(e~^,«)* + A(e^,t)*]) r^} , (17) 



X {Rc(A(Z)/,t/50:f [A(gi?,«)* - A(g^,t)*]) 

+ Im(A(Z)/,f /50:f[A(g«,«)* + A(g«,t)*]) r^} , 



X {Rei[A{eL,t)A{eR,ur - A(e^, x.)A(e«, i)*]/,t/eVet/e1) 

+ Im([A(e^, t)A{eR, u)* + A(e^, u)AieR, ty]f^:fi^f^*f^) r^} , (19) 

with — 1 < Vt < 1 [(^r)^ + (^L )^ ^ 1] being the degree of transverse beam po- 
larizations. The selectron propagators are A{eL,R,t) — i/{t — ?7t.|^^) and A{eL,R,u) — 
i/{u — ^), with the Mandelstam variables t = (p^. — Pe-Y and u — {p^^ — Pe-Y- For 
a centcr-of mass energy ^/s far beyond the Z°-threshold, the Z°-width can be neglected 
in the propagator A(Z°) = i/{s — m\). The kinematical factors in Eqs. (17)-(19) are 

r\ = m^.{[(i_ • s^)(i+ • p^J + (i_ • p^XU ■ s')\{Pe- ■ Pe+) 

+ [{pe- ■ S'){pe+ -pj + (pe- ' Pxi){Pe+ ' s") - (Pe" ' Pe+){Pxi ' " *+)} , (20) 



A = £^.,.m;,Ji>^p^+s^-(i_-p^J+i^ip^p^+p^^(i+.s^) 

+t^-t\pl^s'''^{p,- ■ p^J + e^t^f^.pl^ip,^ ■ s')] , (21) 

where £0123 = —1- The polarization vectors t± are given in Eq. (A. 2), and the neutralino 
polarization vectors are given in Eq. (A.5), for details see Appendix A. 

In [34], the properties of Epj, (16) have been discussed in detail. In contrast to 
unpolarized and longitudinally polarized beams, it contains no contributions from pure 
selectron exchange, Sp(eLei,)T and Sp(ei?eR)T- However, transverse beam polarization 
leads to e^-e^ interference Sp(eiej:j)T, which is absent for unpolarized and longitudinally 
polarized beams [12]. In addition, there is no contribution Ep(ZZ)T owing to the Ma- 
jorana character of the neutralinos. Note that in the high energy limit rrie/ ^/s —>■ 0, the 
neutralino polarization is proportional to the product of transverse beam polarizations 

Xip OC ' 
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The contributions to the amphtude squared, Eq. (12), from the neutrahno decay are 



cos^ Ow 



{L) + R)){2\0:^^\\p^.pj)[m. 



2 



(22) 



■'D 



COS- Ow 

"L\2 
nj 



4){io:?r-x. 



Xj 
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b,fi 1/ , 



[Re(0;^)2 - lm(0;^)>,, 



+ ml) {Pfs\. )-{pz- . ) (Px.- -P/)] } ■ 



(23) 



Different ways to calculate the squared amplitude |Tp for the decay chain x° — > x°^°, 
followed by Z° — > //, with complete spin correlations in the spin density matrix for- 
malism, have been given in detail in [29,33,43]. In contrast to previous calculations, we 
directly obtain the terms D and for the entire neutralino decay chain by summing 
over the internal polarization vectors of the propagating boson in the narrow width 
approximation. 



4 CP asymmetries 

For the combined process of neutralino production and decay (l)-(3), the spin correla- 
tions Sp in the amplitude squared |Tp (12), allow us to define several CP sensitive 
observables. These are asymmetries in the azimuthal angular distribution of the final 
state fermions from the Z° boson decay (3). The CP asymmetries are based on different 
triple product correlations, which we classify in the following according to the beam po- 
larization. We will define CP asymmetries for transversely as well as for longitudinally 
and/or unpolarized beams. 

4.1 Transverse beam polarizations 

We define a CP observable for transverse beam polarizations which is based on the T-odd 
correlation 

Ot = t+ ■ ipf X pe-){t- ■ Pf) + L ■ ipf X Pe-){t+ ■ Pf) = sin(77 - 2(f) f) , (24) 

where p^- is the unit vector of the e" beam, pf the unit vector of the final fermion / in 
Z° — > //, 0/ is the azimuthal angle of /, and the constant = 0+ -|- 0_ is the sum of 

the two azimuthal angles of the polarization vectors of the positron beam, t+, and the 
electron beam, t_. The equality Ot = sm{r] — 2(f) f) in Eq. (24) follows from the specific 
parametrization of the momenta in the center-of-mass system, for their definition see 
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Appendix A. Note that the T-odd correlation Ot is contained in the kinematical factor 
Eq. (21), entering in the neutrahno polarization terms Sp^., given in Eqs. (17)-(19). 
The CP asymmetry is then defined by 

T ^ N[Ot > 0] - N[Or < 0] 
^ N[Ot > 0] + N[Ot < 0] ' ^ ^ 

where the number of events with Ot > 0{< 0) is denoted by N[Ot > (< 0) ]. With the 
definition of the cross section, Eq. (14), C da = dN, where C is the integrated luminosity, 
and the amplitude squared (12) with the decompositions of the quantities P and Ep, 
Eqs. (15) and (16), wc obtain for the CP asymmetry 

T ^ /sgn[0T] iTprfLips ^ /sgn[0r] S^,t dUps' 

f /|T|2rfLips /(P,,poi + Pi) D rfLips' ' ^ ' 

summed over b = 1,2,3, and with rfLips' = dLips/{ds^^dsz) (B.l), where we have al- 
ready used the narrow width approximation for the propagators (B.5). In the numerator 
of Aj, Eq. (26), the spin correlations Spj.S5) are only non- vanishing if weighted with 
sgn[C'T], as they include Or — sm{r) — 2(f)f). In the denominator of Eq. (26), we have 
used J Pt D dLips' = since the cross section is independent of the transverse beam po- 
larizations. Also the terms of the amplitude squared which depend on the polarization 
vanish for an integration over the entire phase space. The asymmetry is sensitive 
to CP violating couplings in the production only, which enter Sp j,, see Eq. (16). That 
means Aj — for the production of a pair of equal neutralinos e'^e" — > x°x°. In addition, 
the asymmetry Aj is proportional to the product of transverse beam polarizations 

A] oc T,''pT oc VtVt > (27) 

and thus vanishes if one beam is unpolarized. 

For the measurement of Aj — ~Aj, the charges and the fiavors of the final fermions 
/ and / have to be distinguished. For f — e, n this will be possible on an event by 
event basis. For / = r it will be possible after taking into account corrections due to 
the reconstruction of the r momentum. For f = q the distinction of the quark flavors 
should be possible by flavor tagging of the heavy quarks q = b. c [46,48]. The distinction 
of the heavy quark charges can be accomplished with very good precision in the case 
of semi-leptonic decays of the hadrons. For the majority of b and c-jets, it also can be 
accomplished by the reconstruction of the vertex charge in the cases where the hadrons 
decay non-leptonically [45,48,49]. However, b and c tagging will be essential. 

The asymmetry for flnal quarks is always larger than that for final leptons Aj, 
due to the dependence of Aj on the Z'^ff couplings [29,33,43] 

AT ^ ~ ^ AT _ + ^bjc) ^bjc) .T /^ox 

^;,) ^ 6.3(4.5)x^J, (29) 
which follows from Eqs. (22), (23) and (25). 
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4.2 Longitudinal or unpolarized beams 



For longitudinal or unpolarized beams, the triple product [36] 

CL=Pe--(p7xP/) , (30) 

can be used to define the CP asymmetry 

L _ /sgn[0L] |Tp dLips _ JsgnjOL] + S^p,J ciLips^ 

/|T|2rfLips /(P,„poi + Pl) c^Lips' ■ ^ ^ 

This CP asymmetry has been analyzed in detail in [33]. Due to the different angular 
dependence of the triple product Ol, Eq. (30), and the CP-odd terms in Ep.^, we have 
/ sgn[Oi;,] Epy dLips' = 0. This means that the asymmetry does not depend on 
the degree of transverse beam polarization. In contrast to the asymmetry Aj, Eq. (26), the 
asymmetry A^ is also sensitive to the CP violating couplings in the neutralino decay Xj 
X°Z° [33]. Due to the dependence of on the Z'^ff couplings, we again have [29,33,43] 

Al^,) ~ 6.3 (4.5) x^^ (32) 



4.3 Measurability of the CP asymmetries 

For a measurement of non- vanishing phases in the neutralino sector at the ILC, not only 
large CP asymmetries but also large cross sections are required. Since the cross section 
appears in the denominator of the CP asymmetries Aj, Eq. (26), and A^, Eq. (31), we 
are often confronted with a situation where large cross sections lead to small asymmetries 
and vice versa. In order to estimate whether the CP asymmetries can be measured at the 
ILC, we consider their statistical significances [29, 33, 43] . 

The significance for the CP asymmetry Aj with transversely polarized beams is defined 

by 

Sj = \A]\^J/: aie+e- ^ x?X?)unpoi BR(x° - X°^°) BR{Z^ ^ //) . (33) 

Note that only the unpolarized cross section a{e'^e~ — > XiXj')unpoi appears in Eq. (33), 
since the cross section is independent of transversal degrees of beam polarizations. We 
define the statistical significances for the asymmetries A^ by 

= lA'fl^C a{e+e- ^ x°X°) BR(x° - x^^o) BR{Z^ ^ //) . (34) 

Here, the cross section depends on the degree of longitudinal electron and positron beam 
polarizations. Note that Sj and Sj, in short Sf, are larger for f = b,c than for f = £ — 
e,/i,r with [29,33,43] 

Sb ~ 7.7 X Si and S^ ~ 4.9 x S^ , (35) 

caused by larger asymmetries (28), (33), and by larger branching ratios BR(Z° — > bb) ~ 
1.5 X E£BR(ZO ^ ii), and BR{Z^ cc) ~ 1.2 x E^BR(ZO ^ ££), £ - e,/i,T, with 
E^BR(ZO ^££) -0.1 [20]. 
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For an ideal detector, a significance of, e.g., Sf = 1 implies that Af can be measured 
at the statistical 68% confidence level. However, the definition of our theoretical signifi- 
cance is based on statistics only, and does not include detector and particle reconstruction 
efficiencies. Our significances are thus upper bounds to judge on the feasibility to measure 
the CP asymmetries. In order to predict the absolute values of confidence levels, detailed 
Monte Carlo analyses including detector and background simulations with particle iden- 
tification and reconstruction efficiencies would be required. A Monte Carlo analysis for 
a CP asymmetry in the production and decay of neutralinos with longitudinal polarized 
beams has been carried out in [35]. However, such an analysis is beyond the scope of the 
present work. We only estimate how the detection rates for b- and c-quark jets have to 
be modified. Using vertex detectors, fiavor tagging of b- and c-quarks is possible and the 
corresponding efficiencies and purities have been studied [45] . It has been shown [46] , that 
6-quarks (c-quarks) can be identified with an efficiency of 50% (50%) at a purity of 90% 
(80%) for e+e" qq at ^/s = 500 GeV. This would lead to a reduction of our statistical 
significance Si, (Sc) by a factor of 0.64 (0.57). 



5 Numerical results 

We present numerical results for the CP asymmetry Aj, Eq. (25), with transverse beam 
polarizations, and for the CP asymmetry A^, Eq. (31), with longitudinal beam polar- 
izations. For neutralino production e'^e~ — > XiX2(3) decay X2{3) ~^ Xi-^^? ~^ 
we study the dependences of A^ and the cross section on the gaugino mass parameter 
Ml = iMile**^^! and on the higgsino mass parameter fi = l^le**^^, as well as on the beam 
polarizations V{e~) and V{e'^). Throughout this study we take ^/s = 500 GeV. The 
values of the CP asymmetries for final leptons Aj, and those for final c-quarks Aj, are 
related to by Eq. (28). Finally we compute cross sections and asymmetries which are 
accessible with longitudinally polarized beams, and compare them with the results for 
transversely polarized beams. 

For the calculation of the neutralino widths and branching ratios, we include the 
two-body decays [29] 

X° ^ ej, i^ei^e, x'lZ', xfW^, , (36) 

with n = R, L for i = e, /x, and n = 1, 2 for £ = r. The Higgs mass parameter is chosen 
niA — 1 TeV. With such a high value for ttia explicit CP violation in the Higgs sector is 
not important for the lightest Higgs state [47]. In the stau sector, we fix the trilinear 
scalar coupling parameter = 250 GeV. To reduce the number of free parameters, we 
take \Mi\ = (5/3) M2tan^6'iy, inspired by gaugino mass unification. 
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(a) a(e+e- ^ [fb] (b) [%] 
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VT[e~) Trie-) 

Figure 2: Neutralino production e'^e" X1X3 decay Xs ~^ ~^ bb, at 

y/s = 500 GeV for Scenario A given in Table 1. Phase dependence of (a) the cross section 
^(e+e- ^ = ^{e+e- ^ x'lXl) ■ BR(x^ ^ " BR(ZO ^ 66), and (b) the CP 

asymmetry Aj for transverse beam polarizations PT(e~) = 0.9, VT{e~^) = 0.6. Contour 
lines in the VT{e~)-VT{e~^) plane for (c) the CP asymmetry Aj, and (d) the statistical 
significance 5"^, with 0Mi = 0m ~ 0' ^^"^ integrated luminosity C = 500 fb^^. 
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Scenario 


M2 




tan P 




men 






X3 




A 


380 


120 


50 


250 


210 


107 


119 


205 


114 


B 


250 


200 


7 


240 


220 


119 


163 


213 


161 


C 


280 


340 


3 


350 


250 


141 


236 


349 


234 



Table 1: Input parameters M2, tan/3, mg^, and mg^ for Scenarios A, B, and C, 
with I Ml I = (5/3) M2ta,n'^ 9w The neutrahno and chargino masses are calculated with 
4>Mi — X 4>ix — ^- All mass parameters are given in GeV. 



5.1 Neutralino X1X3 production and decay Xs ~^ Xi^^ ~^ Xi^^ 

In Fig. 2a, wc show contour lines of the cross section a{e~^e~ — >■ XiXibb) = a{e~^e~ 
X1X3) ■ BR(x3 Xi^^) ■ BR(Z° bb) in the (p^r'pMx plane for Scenario A, see Table 1. 
The production cross section a{e^e^ XiXlDunpoi as well as the decay branching ratio 
BR(x3 — > Xi^^) (which varies between 5% and 15%) strongly depend on (pMi but only 
mildly on 0^, and their minimum values are obtained for — tt and (p/^ — 0. In Fig. 2b 
we show the phase dependence of the CP asymmetry Aj, Eq. (25). The asymmetry Aj 
reaches its maximum of about 36% at (pMi = 1.057r and 0^ = I.Stt. The maximum of the 
asymmetry Aj is close to = ^- The reason for this is that the cross section, which 
is the denominator of the asymmetry, has a minimum there. We remark that the CP 
asymmetry is asymmetric with respect to the transformations — > 27r — and 4>^ — > 
27r— 0^, while any CP-even observable, e.g. the cross section in Fig. 2a, must be symmetric 
under these transformations. Fig. 2c shows the dependence of Aj oc Vt{g^) ■ Vt{g^) on 
the degrees of transverse polarizations of the electron and positron beams for Scenario 
A with = ^ ~ ^- '^^^ maximum value of the asymmetry .4.^ ~ ±40% is 

attained for the maximal degree of polarizations Vrie.^) — ±1 and Vt{g~) = ^1- In 
Fig. 2d we show the statistical significance S},, Eq. (33), for an integrated luminosity of 
C = 500 fb~^. A measurement of the asymmetry thus requires a high degree of the beam 
polarizations. In the following we assume that a polarization of 90% for the electron 
beam, and 60% for the positron beam is feasible, and study the |/x|-M2 dependence of 
the asymmetry Aj. 

In Fig. 3a, we show contour lines of the production cross section a{e'^e'~ X1X3) in 
the \ijL\-M2 plane, with = x — 0- ^he other MSSM parameters are fixed as 

defined in Scenario A of Table 1. In Fig. 3b we show the contour lines of the branching 
ratio BRixl ^ X?^°)- It is as large as 90% for M2 ~ ~ 170 GeV. For |^| > M2 
and 1/^1 > 200 GeV the decay xt ~^ Xi'^'^ dominates. Fig. 3c shows the corresponding 
CP asymmetry A^ in the \iJi\-M2 plane. The asymmetry A^ attains large values in the 
parameter region where the lighter neutralino states x° 2 higgsino-like. For instance, 
for \ix\ ^ 140 GeV and M2 ^ 420 GeV, the asymmetry is .4^ —20%. In addition, 
we show in Fig. 3d the cross section (j[e~^e~ XiXibb) (= o"(e"'"e~ X1X3) ' BR(x3 — >■ 
X?^°) •BR(Z° ^ bb)). For ^ 150 GeV and M2 ^ 170 GeV, we obtain the largest values 
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Figure 3: Contour lines in the \fi\-M2 plane of (a) the production cross section a{e~^e~ — > 
X1X3)) (b) the branching ratio BR(x3 — > XiZ)^ (c) the CP asymmetry Aj, and (d) the 
total cross section a{e~^e~ XiXi^b), for Scenario A given in Table 1, with (pMi = X' 
(j)^ = 0, at a/s = 500 GeV with transverse beam polarizations Vrie") = 0.9, Vrie'^) = 
0.6. In the gray shaded area m^± < 104 GeV. The dashed line in (a) indicates the 
kinematical limit m^o + m^o = ^/s, and the dashed-dotted line in (b) indicates the limit 
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up to 7 fb, where the asymmetry Aj reaches 8.5%. With these values of and a, the 
statistical significance, Eq. (33), is S*^ ~ 5, for an integrated luminosity of C — 500 fb"-*^. 

Now we discuss in neutralino production e~^e~ — >• X1X3 ^^'^ decay Xs ~^ Xi^°) 
— > bb, for Scenario B, see Table 1, with a small value for tan/3 — 7. In this sce- 
nario, the neutralino Xs is a strong mixture of gaugino and higgsino components, and 
BR(x3 Xi^^) = 100%. In Fig. 4a we show contour lines of the production cross section 
a{e^e~ X1X3) the 0^-0Mi plane. The production cross section ranges from 25 fb, 
for 4>Mi = TT and (p^ = 0, to 54 fb, for =0/^ = 0. In Fig. 4b, we show the total cross 
section a{e'^e~ — > XiXibb), which varies between 4 fb and 8 fb. Contour lines for the CP 
asymmetry Aj are shown in Fig. 4c. As can be seen, Aj can attain large values even for 
CP phases close to the CP conserving case. For instance, for (J)^ = and (pMi = O.Stt, 
the asymmetry reaches about —10%. The corresponding statistical significance, shown in 
Fig. 4d, reaches = 5, for £ = 500 fb^^ 

5.1.1 Comparison of transversely and longitudinally polarized beams 

In this subsection we compare the asymmetry Aj, Eq. (25), that is obtained with trans- 
verse beam polarization with the asymmetry A]^, Eq. (31), for longitudinally or unpolar- 
ized beams. Fig. 5 shows the contours of the asymmetry A]^ and the total cross section 
(j(e+e~ — > x?x?6fc) for longitudinally polarized beams in the |//|-M2 plane for scenario A, 
see Table 1. The asymmetry A^ can reach ±10%, while the total cross section can reach 
15 fb. By comparing these results for longitudinal beam polarizations with the results of 
Fig. 3 for transverse beam polarizations, we find that in general for higgsino-like scenarios 
the statistical significance of asymmetry Aj is larger than that for asymmetry Aj^. The 
reason is that for higgsino-like scenarios the leading contribution to the asymmetry Aj is 
due to the e^-CR interference term, Eq. (19), which is absent for longitudinally polarized 
beams. Clearly, for larger values of the selectron masses m^^ and m^^, the significance 
of the asymmetry A^ is reduced due to a suppression of the e^-Cij interference term, 
as these masses enter in the propagators, see Fig. 1. On the other hand, for gaugino- 
like scenarios, the statistical significance of asymmetry A^, is larger than that for Aj, 
since then the cr-cr and the cl-gl terms, which arc only present for unpolarized and 
longitudinally polarized beams, give the dominant contribution to A^. 

In Fig. 6, we show the phase dependence of the production cross section a{e^e~ — > 
X1X3)) the total cross section ^(e+e^ XiXi^^)i the asymmetry A}j^ and the corresponding 
significance for longitudinally polarized beams {V^-,V^+) = (+0.9, —0.6) in scenario A of 
Table 1. The asymmetry A^ can reach ±15%, the total cross section can be as large as 
(7(e+e~ — > XiXibb) — 16 fb, such that the significance can reach S^' — 12. In comparing 
these results. Fig. 6, with the results for transversely polarized beams, see Fig. 4, we find 
that the asymmetry A]^ can be measured with higher statistics compared to Aj. The 
reason is that the high degree of longitudinally polarized beams, {V^-,V^+) = (0.9, —0.6), 
enhances both the production cross section and the asymmetry, whereas transversely 
polarized beams do not change the cross section, see Fig. 4. 
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Figure 4: Neutralino production e^e~ — X1X3 and decay Xi^^j ~* at 

y/s = 500 GeV with transverse beam polarizations Trie') = 0.9, Vt^^^) = 0.6, for 
Scenario B, given in Table 1. Contour lines in the (t)pr<pMi plane for (a) the production 
cross section a{e'^e~ X1X3)) (b) the total cross section a{e'^e~ XiXibb), (c) the 
CP asymmetry Aj, and (d) the statistical significance S^, for an integrated luminosity 
£ = 500 fb-^ 
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Figure 5: Neutralino production e+e~ X1X3 and decay Xs ~^ ~^ bb at 

■s/s = 500 GeV with longitudinal beam polarizations Vii^') = 0.9 and Viie^) = —0.6. 
Contour lines in the |yu|-M2 plane for (a) the CP asymmetry A]^, and (b) the total 
cross section (T(e^e~ XiXibb), for Scenario A, given in Table 1. The branching ratio 
BR(x3 -> XiZ°) is shown in Fig. 3(b). 

5.2 Neutralino X1X2 production and decay X2 ~^ Xi^^ Xi^b 

In Fig. 7, we show the phase dependence of the CP asymmetry Aj, Eq. (25), and the 
cross section for neutralino production e~^e~ — ^ X1X25 decay X2 ~^ x\Z^ 1 ~^ bb, 
for Scenario C, see Table 1. The asymmetry reaches a maximum value of about 4.5% 
for (</)Mi,0Ai) ~ ("~f)f)- The production cross section a{e^e~ X1X2) varies from 
19 fb for (0Mi,0^) = (0,0), to 72 fb for {4>Ml,(f^^i) = (7r,vr). The branching ratio is 
maximal BR(x^ x?^°) = 100% for (p^ < OAn and 0^ > 1.67r. For OAtt > (p^, < I.Gti, 
the X2 decay channels X2 ~^ ^R^ kinematical accessible, which leads to a strong 0^ 
dependence of the cross section cT(e+e~ XiXV^) around the threshold region. In the 
region where BR(x2 Xi-^°) = 100% (0^ < 0.47r and 0^ > 1.67r) the total cross section 
o"(e^e~ — » XiXi^^) varies from about 4 fb to 10 fb, see Fig. 7b. For Scenario C, with 
phases (0Mi, 0^) = (— f , f ), the significance is Sb ~ 2, for C = 500 fb"-*^. 

In Fig. 8, we show the asymmetry A^ and the cross section for longitudinally polarized 
beams {V^-,V^+) = (—0.9, +0.6) in scenario C of Table 1. The asymmetry A]^ can reach 
±10%, and the total cross section can reach a{e^e~ XiXi^b) = 25 fb. Again the strong 
0^ dependence of cr(e~'"e~ XiXi^b) in the regions where 0^ ~ 0.47r, 1.67r is because 
the decay channels X2 ~^ ^R^ kinematical accessible beyond this threshold. The 
significance to measure for (0Mi, 0^4) = (~f ? f ) is S'^ = 7 (£ = 500 fb~^). Once more 
we find that for longitudinally polarized beams, the asymmetry A}^ can be measured with 
higher statistics than the asymmetry . 
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Figure 6: Neutralino production e+e~ X1X3 and decay Xs ~^ Xi^^j ~^ bb at 
y/s = 500 GeV with longitudinal beam polarizations ^^(e") = 0.9, Viie'^) = —0.6, for 
Scenario B, given in Table 1. Contour lines in the (t)pr<pMi plane for (a) the production 
cross section a{e'^e~ X1X3)) (b) the total cross section (T(e+e~ — >• x\x\hb), (c) the 
CP asymmetry and (d) the statistical significance S"^, for an integrated luminosity 
C = 500 fb-^ 
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Figure 7: Neutralino production e^e~ X1X2 decay X2 ~^ Xi^^i ~^ ^ 
y/s = 500 GeV with transverse beam polarizations Pr(e~) = 0.9 and Pr(e^) = 0.6. 
Contour lines in the (p^-4>Mi plane for (a) the CP asymmetry A^, and (b) the total cross 
section cr(e+e~ XiXibb), for Scenario C, given in Table 1. 

Finally, we have scanned the MSSM paramter space, and have found that in general 
the CP asymmetry Aj is smaller for X1X2 production than for X1X3 production. Note that 
in the parameter region where M2 > for which the largest values of Aj are obtained 
for X1X3 production, the two-body decay X2 ~* Xi^^ is kinematically forbidden. 

6 Summary and conclusions 

We have studied the impact of the CP violating MSSM phases 0Mi and 0^ on neutrahno 
production e~^e~ XaX^ decay x° Xn-^°) ^ //• transversely polarized 
beams, we have defined the CP observable Aj, which is an asymmetry in the azimuthal 
distribution of the final fermions, and is based on the triple product correlation Eq. (24). 
The CP-asymmetry depends bilinearly on the beam polarizations Aj oc V^Vj^, and 
receives CP sensitive contributions from spin-correlations in the neutralino production 
process. Due to the left-right structure of the boson couplings to the final fermions, 
the asymmtry is largest Aj = 6.3 x Af for the hadronic decay —>■ bb. 

In a numerical analysis for X1X2 X1X3 production, we have shown that the asym- 
metry can be as large as Aj = 30%, for Trie') = 0.9 and Vrie^) = 0.6. The significance 
can be as large as = 5, with an integrated luminosity C = 500 fb~^. We have compared 
the cross sections and asymmetries with those which are accessible for longitudinally po- 
larized beams. For X1X2 production, we have found larger significances for the asymmetry 
A^ with longitudinally polarized beams. Also for X1X3 production the asymmetry Ah is 
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Figure 8: Neutralino production e^e~ X1X2 decay X2 ~^ ~* bb at 

■s/s = 500 GeV with longitudinal beam polarizations PL(e~) = —0.9 and PlIc"*") = 0.6. 
Contour lines in the 0^-0 a/i plane for (a) the CP asymmetry A};, and (b) the total cross 
section a{e~^e~ — > XiXi^b), for Scenario C, given in Table 1. 

in general accessible with a larger statistical significance compared to the asymmetry Aj 
based on transverse beam polarizations. This is mainly because longitudinal beam polar- 
izations can greatly enhance the production cross section. There are, however, parameter 
regions, e.g. for M2 > 300 GeV and < 200 GeV, where the asymmetry Aj in X1X3 
production can be up to about a factor five larger than the asymmetry A^ 

We conclude that for neutralino X1X2 production and subsequent two-body decay into 
the Z° boson, longitudinally polarized beams give more statistics than transversely po- 
larized beams {S^ > Sj) to study the corresponding CP asymmetry in the neutralino 
sector of the MSSM. However, for X1X3 production, there are parameter regions where 
transversely polarized beams are advantageous since Sj > S^. We emphasize that the 
asymmetries for longitudinally und transversely polarized beams are independent of each 
other. They rely on different triple product correlations, and depend on different interfer- 
ence channels in the neutralino production process. Therefore both options of polarized 
beams should be considered to determine the phases 4>Mi and 0^ by measurements of the 
asymmetries Aj and A^ at the ILC. 
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Appendix 

A Momentum and polarization vectors 

We introduce a coordinate frame in the laboratory system by choosing the 2;-axis along 
the momentum vector of the e~ beam, and x and y according to a right-handed coordinate 
system: 

=^6(1,0,0,1) and p^+ = £;ft(l,0,0,-l) , (A.l) 

with the beam energy Ei, — -\/s/2. Then the 4- vectors of the transverse polarization of 
the electron and positron beams are in the x-y plane 

t± = (0, cos 0±, sin 0±, 0) , (A.2) 

where the azimuthal angles 0+ and 0_ describe the orientation of the transverse beam 
polarizations, which are fixed at any value G [0, 27r). The four-momenta of the 

neutralinos are 

Pv = l{ — ^, COS sin 6', sin sin 6', COS 6') , 
X.3 q 

E 

Pxi ~ Qi^^, — cos(j)sm9, — sm(j)sm9, — cos9) , (A. 3) 

where 9 is the scattering angle and (f) the azimuthal angle of the production process. The 
energies and momenta of the neutralinos are 



^ ^ x^o)^ ^ ^ (A.4) 

2^5 2^5 



with A(a,6,c) = ^-h"^ + — 2{ah + ac + hc). The three spin-basis vectors s^'^^, h = 1,2,3, 
of neutralino x° are 

0,^1 ^ = (0,-cos(/)cose,-sin(/)cose,sin^) , 



^x. 



s^'^ = ^— I g, — jl. ) = ^^(— ,cos0sin6',sin0sin6',cos6l) , (A.5) 

m^. V g V "^x.- ^xi 
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with s^^ ■ Py^. = and {s'^.s^^.) = —Sab, such that {^xj^ ^xj^ ^Xj^ build a right-handed 
system. 

For the two-body decay x° — > x°^, the decay angle Oil{pxj,'Pz) is constrained by 
g-j^^max _ qO^q fQj. ^ where 5° = A^(m^^., m|, m^^)/(2m2) is the neutralino momen- 
tum if the Z° boson is produced at rest. In this case there are two solutions [29,33,43] 



{ml. + ml- ml Jq cos 01 ±E^.^X{ml.,m%, ml J -iq'^ ml {1-cos'^ei) 
'^^ ' " 2g2(i-cos2 ^i) + 2m2. ■ ^ '^^ 

If g° > q, the decay angle 9i is not constrained and there is only the physical solution 
\Pz\- The momenta in the laboratory system are 

"i| = IpJ^I sin^i COS01, sin^i sin^i, \p^\cos6i) , (A. 7) 

= IpjI sin 62 cos 02, sin 82 sin 02, cos 62) , (A. 8) 



Pf 



ml , , 

^ (A.9) 



2(E|-b-'±|cosM ' 
with 92/-iPxj,Pf) decay angle OD/-{pz-iPf) given by 

cos 9d — cos ^1 cos 82 + sin ^1 sin 82 cos(02 — 0i) ■ (A. 10) 

B Phase space 

The Lorentz invariant phase space element for the neutralino production (1) and the decay 
chain (2)-(3) can be decomposed into the two-body phase space elements [29,33,43] 

dLips(s,p^.,p^„,P/,p/) = 
1 



(27r)^ 



dLips(s,p^,,P;^.) ds^^ ^dLips(sx.,p^„,p|) dsz dUps{sz,Pf,Pf) , (B.l) 



dUps{s,Px„Pxj) = 4(2^)2^ ' (^-2) 

1 Ip^l"^ 

dUps{sxj,Pxn,Pz) = T^TT^ 01 F± F |->±|| ^^1 ^ (B-3) 



2(27r)2 2|E| gcos^i-S^, |p# 
. dO.^ 

2(27r)2 m| 



dUps{sz,Pf,Pf) = ——7^ o?n2 , (B.4) 



with s^. = pI^, sz = pI and dfli — sin^j d9i d(j)i. We use the narrow width approximation 
for the propagators 

/ ds,, I \A{Z')\'dsz = ^ . (B.5) 
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The approximation is justified for T^./m-^. -C 1, wliicli liolds in our case witli T^^. < 
O(lGeV). Note, liowever, tfiat tfie naive C(r/m)-expectation of the error can easily 
receive large off-shell corrections of an order of magnitude and more, in particular at 
threshold or due to interferences with other resonant or non-resonant processes. For a 
recent discussion of these issues, see [53, 54] 
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